Alternative mechanism of avoiding the big rip or httle rip for a scalar phantom field 
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Depending on the choice of its potential, the scalar phantom field (j) (the equation of state pa- 
rameter w < —1) leads to various catastrophic fates of the universe including big rip, little rip and 
other future singularity. For example, big rip results from the evolution of the phantom field with 
an exponential potential and little rip stems from a quadratic potential in general relativity (GR). 
By choosing the same potential as in GR, we suggest a new mechanism to avoid these unexpected 
fates (big and little rip) in the inverse- -R gravity. As a pedagogical illustration, we give an exact 
solution where phantom field leads to a power-law evolution of the scale factor in an exponential 
type potential. We also find the sufficient condition for a universe in which the equation of state 
parameter crosses w = —1 divide. The phantom field with diff'erent potentials, including quadratic, 
cubic, quantic, exponential and logarithmic potentials are studied via numerical calculation in the 
inverse- i? gravity with correction. The singularity is avoidable under all these potentials. Hence, 
we conclude that the avoidance of big or little rip is hardly dependent on special potential. 



I. INTRODUCTION 



<«J 'i It has been widely accepted that astrophysical observations might favor a dark energy component with " superneg- 
' ative" equation of state (EOS) parameter w = p/q < —1, dubbed as phantom or super quintessence, where p is the 
f^, pressure and p is the energy density [ll, [3l ■ The phantom model with Born-Infeld type Lagrangian has been proposed 
Q and its generalization to 0{n) symmetry has been done in Ref.Q. The phantom component leads to a catastropic 
big rip singularity that is characterized by divergences in the scale factor a, the Hubble parameter H Q and its 
J> ■ time-derivative [gl at the finite future t = It has been studied that various types of singularity might occur in 
•/"J , phantom scenarios [3, |1| . Barrow and his collaborators 0, [13] have investigated general characterization of sudden 
■ cosmological singularities and the classical stability of cosmological solutions containing these singularities. Some 
, mechanisms have been proposed to avoid these catastropic singularities. For example, it has been shown that these 
singularities in the future of the cosmic evolution might be avoided if we take suitable potential term [lT| - [l3| , quantum 
effects fT3] and modified gravity ^5j. If one considers an additional interaction between the phantom field and the 
background, the big rip singularity can be avoided in the late-time evolution of the universe [l6j . In brane- world 
models, a Gauss-Bonnet term is provided for the bulk, whereas phantom field is present on the brane. If the dark 
energy is modelled by a phantom-generalized chaplygin gas, the evolution of universe will avoid the big rip [T7l - [l9j . 
In the torsion cosmology, the scalar mode of torsion could be considered as a phantom field which will derive the 
universe in an oscillating fashion with an accelerating expansion at late-time [2Cll424| . However, w < —1 as the scale 
5^ ■ factor a{t) extends into the future is necessary condition [3,3, but not sufficient one for a future singularity. Recently, 



X 

d I the little rip scenario has been proposed in [25[, in which w < —1 but u> ^ — 1 asymptotically, such that there is no 
future singularity. Such model will nevertheless lead to a dissolution of bound structures at some time in the future. 
Especially, the viscous fluid can produce a little rip cosmology as a purely viscosity effects [26| . 

On the other hand, there is an interesting way to give rise to acceleration: a modification to the Einstein-Hilbert 
action involving new terms of inverse powers of the curvature scale (in the following text we call it inverse- i? gravity), 
of the form y/—gR^" (with n > 0) |27|. In general f{R) gravity (even without scalar field), the divide w = — 1 
crossing behavior is shown in Refs. 



31j . Recently, Du et al. [32[ discussed the behavior of quintessence field in 
the inverse- i? gravity and found that the evolution of w with cosmic time can cross the divide w = —1 and w will be 
less than —1. In this paper, we put forward a new mechanism to avoid the catastropic big or little rip by means of 
the investigation of phantom field in the inverse- i? gravity. It is worth noting that crossing w = —1 divide happens 
approximately in a few e-fold expansion when the inverse- i? term is still sufficiently large and the decrease of the 
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energy density resulting from the crossing of w will make the catastrophic fate avoided. 

This paper is organized as follows. In Sec. 2, we give a brief review of inverse- i? gravity and its relation to the 
phantom field in Einstein's general relativity (GR). An exact solution is shown in Sec. 3. In Sec. 4, we investigate 
the mechanism for a single phantom field to avoid big or little rip by an analytical method. In Sec. 5, we present the 
numerical results to avoid big or little rip of phantom depending on different potentials in frame of inverse- i? gravity. 
Finally, we give a brief summary and a further discussion in Sec. 6. 



II. PHANTOM FIELD IN INVERSE- i? GRAVITY 



There are two types of modification of Einstein-Hilbert action: one is ultraviolet modification which is important 
in the early universe (high curvature region) and the other is infrared modification which plays significant role in the 
late universe (low curvature region). The inverse- i? gravity is one of the current models with infrared modifications. 
The inverse- i? gravity seems unable to pass the solar system observations 33], but this problem has been solved by 
adding a scalar curvature squared term to the inverse- iZ gravity [s^, |35|. In Ref. [sT], the following f{R) theory 
was proposed: 

where n,m,Ai,A2 and /? are constants. This model leads to an acceptable cosmic speed-up and is consistent with 
some of the solar system observations. This general case returns to the inverse- i? gravity when the curvature is low 
and Ai = 0. 

In this paper, we consider the following action 

3="^ j d^x^g{R~^ + ^)+ j d^x^g{C^+C^) (2) 

in the FRW universe with scale factor a. The metric is 

ds^ = -df + a^dY? (3) 

where dS^ is the metric of a 3-dimensional maximally symmetric space. £,„ is the Lagrangian of matter and Ccj, is 
the Lagrangian of phantom field with the potential V{<f)^ 

C^^\g^-d,cj,d,cp-V{4>)- (4) 



The Friedmann equation is given by [31| 



jj2^t^ ^iH,H,H) ^ pm - ^^2 + F(0) 
a2 3m2; 3mli 



where H = d/a is the Hubble parameter and the dot denotes the derivative with respect to t, and 
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^ ^ 12(i7 + 2i/2)3^ ' 

' ^{6H^H - + 2HH). (6) 



/32 

The equation of motion for phantom field is 



• dV 

3H(f)- — = 0. (7) 



In this system the phantom field together with the extra geometric term takes effect of dark energy. In the inverse- ii 
gravity, there is an additional J^-term in the modified Friedmann equation (5). To expound the observed evolving 
EOS of the effective dark energy, one can introduce the concept "equivalent dark energy" [s^. We derive the density 
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of equivalent dark energy by comparing (5) with the standard Friedmann equation in GR. The latter can be written 
as 

H^ + ^ = lT^iPra+Pde)- (8) 

Comparing (8) with (5), one has 

Pde = P4>+T, (9) 

where = — </)^/2 + V{4>). It is worth noting that J^-term makes a key difference from the standard model of GR. In 
the next section, we shall study analytically the phantom dynamics of universe in inverse- i? gravity. 



III. POWER-LAW SOLUTION IN A MODEL WITH EXPONENTIAL TYPE POTENTIAL 

To find exact solutions is an important but difficult topic in such a highly non-linear system (7) and (10). The most 
interesting exact solution is a power-law one because it has the extra advantage that the equations for the generation 
of density perturbations can also be solved exactly in the cosmological inflation [36]. In a spatially flat universe, the 
Friedmann equation (5) reduces to 



3m2, 3™^^ 



during the phantom dominated stage. Note that crossing w — —1 divide happens approximately in a few e-fold 
expansion when the inverse- i? term is still sufficiently large and the i?2-term is still sufficiently small, therefore, the 
/3-term can be neglected in rhs of (6) at the late-time of evolution. We have 
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T(H, H, H) = : ^ (2HH + 15H^H + 2H^ + 6H^). (11) 

In this case, the power-law solution arises when the potential is chosen to take the following form 



l+7exp(-^)l, (12) 



V — Vq exp 
where ^ > 4, 7 and Vq are constants, and 

^ 12(14^3 _ 64^2 _ _ 17)^2 ^_ 4(24^3 _ 64^2.^ 37^ _ 3)^^3^2—12^- 
The equations of motion (7) and (10) then have the particular solution 

a^aot^, (14) 



— =^lnt, (15) 
mpi 



where 



Provided that j > 2 when f > 4, this solution satisfies the condition for accelerating expansion. If we take another 
value of 7 in the potential (12) 



^ 3«4(12-3g+ V3g2-40 

^ 12(14^3 - 64^2 - 70C - 17)^2 + 4(24^3 _ 64^2 _^ 37^ _ 3)^^3^23:4^' 

the system (7) and (10) have the particular solution (14) and (15), but 



Provided that 1 < j < 2 when 4 < ^ < 6, this solution also satisfies the accelerating condition. It is not surprising 
that there is a power-law solution for the dynamical system of phantom because J^-term makes a key difference from 
the standard model of GR. The power-law solution makes a message known to us: the big rip can be avoided in the 
inverse- i? gravity. 
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IV. THE AVOIDING MECHANISM 

Before studying the avoiding mechanism in this dynamical system, we discuss briefly the relation between density 
evolution and EOS. Since the dust matter obeys the continuity equation and the Bianchi identity keeps valid, equivalent 
dark energy itself also satisfies the continuity equation 

^+3H{pde+Peff)^0, (19) 

where Peff denotes the effective pressure of the equivalent dark energy. From (9) and (19), we can express EOS as 

^ ^ ^ ^ _1 _ \^A p^ + ^) (20) 

Pde iH[p^ + F) 

where 

j" = '^"^P'". — \H{2H'^ + H) + 3H(2H'^ - 7HH - H) + 3H(H^ - WH'^H - AH^)] (21) 

and P(p,J-' depend on the rate of expansion. If + < (p^ + T > 0) the equivalent dark energy behaves as 
quintessence (phantom) in GR, i.e., pde decreases (increases) with expansion of universe. 

Next we consider the approach of Ref. 0, where the pressure is expressed as a function of the density in the 
form Peff = -pde - fiPde)- Taking a power law for f{pde), i- e., /(pde) = ^pJe (^ > and 7 > 0),we see that the 
fate of universe will depend on 7 : (i) There is the possibility of little rip but no finite future singularity for 7 < | ; 
(ii) There is a type I (big rip) singularity for i < 7 < 1; (iii) There exists a type III singularity for 7 > 1. Here, the 
singularities are classified according to Ref. [7]. We have 

_ I -noe , / 2' /r,r,\ 

I ^[(2/3-l)A(t,-t)] — , 7>iand7^1 ^ ' 



where h = c = nipi = 1 is taken for the sake of clarity. Substituting (22) into (21), we obtain 

I ~ 2 

7 > i and 7 7^ 1 



) 3M{t)Nit) .1 , , , y^"^) 



2L{t) 



where 



L{t) = (v^+ 2(2/3 -l)(t,-t)[A(2/3-l)(t,-t)]T^)4, (24) 

M{t) = [Ai2l3-l){ts~t)]^, 

N{t) = 2V3(l + 2/3)/3+(2/3-l)(is-t)[^(2/3-l)(t. -i)]T^(46/3-16/32-3) 

+ 4^3(2/3 - ifit, ~ tf[A{2l3 - l)(t, - t)]^ 

+ 4(2/3 - l)3(t, - tf[A{2(3 - l)(i, - t)]^. 

It is easy to prove < ( j" is a negative definite function) for 7 < jcrit, Jcrit = (23 + •\/481)/16, which brings a 
possibility to avoid the big and little rip, and some type III singularity for 7 < jcrit- In the inverse- i? gravity, the 
J^-term makes a key difference from the standard model of GR. 

In the w(tQ) > —I case, w(t) will monotonically increase since is a negative difinite function for 7 < ^crit- In the 
w{to) < —I case, we need an analysis in detail. Combining (7) and (10), we find that 

*=2f^, (25) 
H= a 2 TT • (26) 
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If w{to) < —1 from the current observation, we must demand that EOS parameter cross w = —1 divide for the 
avoidance of big or little rip in the future. The universal argument is that the crossing phenomenon happens at t = tf, 
iff H{tc) — and H{tc) ^ 0. In the phantom case of GR, we have = such that H{tc) = implies H{tc) — 0. 
Therefore, the crossing behavior never happens. In other words, the big or little rip is inescapable in the Einstein's 
GR. In inverse- i? gravity, H{tc) = and H{tc) ^ iff = —2>H(jP and JF ^ —&H(j)(f) at t = tc- Thus, we have proved 
the following theorems. 

Theorem. If w(io) < the avoidance of big or little rip demands J- — —3H(j)^ and ^ —6H(j>(j> at t — tc in the 
inverse- i? gravity where tc is the crossing time. 



V. NUMERICAL CALCULATION FOR THE SCALAR PHANTOM MODELS IN INVERSE-i2 

GRAVITY WITH CORRECTION 

For a minimally coupled phantom field cj) in GR, the necessary conditions have been given for the refrainment of big 
rip or little rip 371 . That is, when V' /V — ?> oo and / VV{V')~^d<l) — > oo, w approaches -1 sufficiently rapidly that big 
rip is refrained. Furthermore, the boundedness of the potential V{(j)) indicates whether or not the model corresponds 
to a little rip expansion. Thus, we examine a specific model with V{4>) — (X) as — >■ oo in more detail. 

In Sec. 4, we have investigated analytically the evolution of universe in the phantom dominated epoch in inverse- 
gravity. For a realistic universe, the pressureless matter is an essential component. However, it is difficult to study 
analytically the evolution of universe when we consider matter component in inverse- i? gravity with correction. 
The behaviors of phantom field with different potentials are investigated by numerical calculation in this section. 
The results show that the big or little rip can be avoided for these models. A spatially fiat FRW universe has been 
assumed in all calculations. 



A. Exponential potential: An example of avoiding big rip 

The form of exponential potential is 

y(0) =/.eexp(^) (27) 
rripi 

where fic is a parameter with dimension of energy density. Let us introduce the four dimensionless parameters for the 
numerical calculations, x = = :p^,b = and cq = Thus, (7) and (10) can be rewritten as 

= fl„^fie-^' - - u{y) + b{x) + c{x), (28) 



where 



x^y" + xx'y' + 3x^y' - v{y) = 0, (29) 

, hUx^x'^ + 2x^x" + lbx^x' + Qx'^) 

= 36(x.- + 2x^)3 ' (3") 

c{x) = 18co(6x^a;' + x^x'"^ + 2x^x"), (31) 

^(y) = (33) 

Here the prime denotes the derivative with respect to the so called e- folds s = Ina and f2m,o = Pm,o/ {^H^ni^i), Hq 
is the current Hubble parameter and Pm,o is the current density of pressureless matter. Generally speaking, cq is a 
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FIG. 1: The evolution of 1 + w,Qde and q with the e-folds s for an exponential potential, where the real, dot, dash lines 
represent 1 + w, fide and q, respectively. 
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FIG. 2: The behaviors of 1 -I- m for the scalar phantom model with an exponential potential in GR (the dash line)and in 
inverse- 7? gravity (the real line). 



little tiny number since the inflation scale /3 is much higher than Hq. If we fix the inflation scale as a GUT energy 
scale of 10" GeV, cq ~ 10-"2_ 

The deceleration parameter can be expressed as 

q = -l-- (34) 
X 

and we introduce the dimensionless density of equivalent dark energy fide for convenience 

= ^J!t 2 - -h^y'^ + <v) + bi^) + (35) 

3-^0%; 6 

Using numerical calculation, we plot w, fide and q as functions of e-folds s for the scalar phantom in inverse- gravity 
in Fig. 1. In Fig. 2, we show the evolutions of w for the scalar phantom model with an exponential potential in 
GR and in inverse- i? gravity. Obviously, the EOS of equivalent dark energy changes from less than —1 to larger 
than —1 such that the big rip singularity is avoidable in the inverse- gravity whereas w is always less than —1 and 
asymptotically tends to a constant which leads to the big rip singularity in GR. It is clear that i?~^-term plays an 
essential role for the avoidance of big rip. The result of numerical calculation is consistent with that of exact power 
law solution Eqs. (14) and (15) for exponential type potential. 
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FIG. 3: The evolution of l + u;, Qde and q with the e-folds s for the quadratic potential, where the real, dot, dash lines represent 
1 + w, fide and q, respectively. The initial values are taken as xo = l,yo = 1, x'q = —1.03, y'o = 0.01. 



B. Quadratic potential: An example of avoiding little rip 

The quadratic potential is widely investigated since any potential around its minimum (if it has a minimum) can 
be treated as the form of quadratic. We take the potential as 

= ^llcp^ (36) 

where /X2 is a parameter with dimension of mass. In this case, (7) and (10) are still reduced to (28) and (29), but 

"(y) = (37) 



v{y) - ^y- (38) 
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In GR, phantom model with quadratic potential may lead to the little rip in the future of cosmic evolution. However, 
we see that the EOS of equivalent dark energy changes from less than —1 to larger than —1 in Fig. 3. It means that 
the little rip is avoidable for phantom field with a quadratic potential in the inverse- i? gravity. 



C. Other potentials 



In this subsection, we study other potentials for phantom field including cubic, quantic and logarithmic potentials. 
In GR, the scalar phantom field with quadratic, cubic and quantic potentials might lead to little rip whereas the big 
rip may occur for exponential and logarithmic potentials [37|. In the numerical calculation, we only need to revise the 
expressions of u{y) and v{y) in (28), (29) and (35) for different potentials, which are shown in Table 1. 



TABLE I: The expressions of u{y) and v(y) for various potentials 
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In Figs. 4-6, we show the numerical results for the cubic, quantic and logarithmic potentials, respectively. We find 
that the results are analogous to the case of exponential or quadratic potential, i.e., the big or little rip is avoidable 
in the inverse- i? gravity. 
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FIG. 4: The evolution of 1 -I- Wjflde and q with the e-folds s for the cubic potential, where the real, dot, dash lines represent 
1 -|- w, Qde and q, respectively. 
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FIG. 5: The evolution oi 1 + w, Qde and q with the e-folds s for the quantic potential, where the real, dot, dash lines represent 
1 -|- KJ, Qde and q, respectively. 



VI. CONCLUSION 

In GR, the phantom field with an exponential or a quadratic potential leads to big or little rip that is characterized 
by the destruction of all bound structures. We show the possibility of avoiding the big or little rip naturally in the 
inverse- i? gravity. For a special exponential type potential, we find an exact power-law solution that makes a message 
known to us: the big rip can be avoided in the inverse- i? gravity. We also find the sufficient condition for a universe in 
which the equation of state parameter crosses w = ~1 divide in the inverse- i? gravity. In GR, the scalar phantom field 
with quadratic, cubic and quantic potentials might lead to little rip whereas the big rip may occur for exponential 
and logarithmic potentials. And then we investigate these potentials for phantom field (j) in the inverse- i? gravity with 
correction. The big or little rip is avoidable under all these potentials. Therefore, we conclude that the avoidance 
of big or little rip is a robust property of inverse- i? gravity with corretion, not rigidly adhered to some special 
potential. 

It is worth noting that crossing w = ~l divide happens approximately in a few e-fold expansion when the inverse- i? 
term is still sufhciently large and the decrease of the energy density resulting from the crossing of w will make the 
catastrophic fate avoided. 

Finally, we give a brief discussion. It is not a shortcoming that we discuss only the background evolution in this 
paper but not consider the cosmological perturbation theory. In fact, Nojiri and Odintsov[3l| have pointed out that 
modified gravity equations of motion are higher-derivative differential equations. The cosmological perturbations in 
modified gravity approaches to those in GR, since the corresponding equations are reduced to second-order differential 
equations. 
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FIG. 6: The evolution of 1 + ui, Qde and q with the e-folds s for the logarithmic potential, where the real, dot, dash lines 
represent 1 + w, fide and q, respectively. 
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